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ABSTRACT 


This thesis reports the results of an experimental investigation into 
the spanwise correlation of the fluctuating longitudinal velocity component 
in the early wake of a flat plate placed parallel to a uniform stream and 
allowed to vibrate about a vertical axis through the leading edge. The 
experimental program provides for the measurement of fluctuating total head 
signals at one fixed and one variable spanwise position on the wake center- 
line one body thickness downstream of the trailing edge. Analyzation of 
these signals gives the spanwise correlation of the wake structure expressed 
as a complex number and the power spectral density of the longitudinal 
velocity fluctuation along the span. In addition the mean local wake 
velocity is measured for each value of the ambient flow velocity. 

Ambient velocities of 8, 11, and 14 fps and a single plate of 4" thick- 
ness and 2" chord, were used for all the tests which were carried out in the 
73" x 9" test section of an open-circuit water Se licis the Hydrodynamics 
Laboratory at MIT. 

An equation is developed which relates the vortex-induced forcing 
moment on the foil to the measured structure of the early wake. Using this 
model, forcing moment coefficients are evaluated as a function of a "self 
excitation" parameter consisting of the ratio of the mean transverse velocity 
of the trailing edge to the free stream velocity. The role of spanwise 
coherence of the wake is examined with respect to its effect on the coupling 


between body vibration and forcing moment. 


Major Conclusion of this study is: 


The increase in the forcing moment coefficient, Me accompanying an 
Oo 
increase in the self excitation parameter is due to increasing correlation 
or coherency of the wake along the span rather than to increasing strength 


of the individual vortices being shed. 
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transverse fluctuating velocity; fps 

mean transverse velocity in the wake; fps 
2BOe "self excitation" parameter; dimensionless 
spanwise fluctuating velocity in the wake; fps 
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cross power density; signal A, signal A 
cross power density; signal B, signal B 


cross power density; signal A, signal B 





TABLE OF CONTENTS 


I INTRODUCTION 


A. Flow Induced Vibrations. . . 


B. Scope of the Investigation . 


II EXPERIMENTAL EQUIPMENT 


Ae Water Tunnel a 


MWmoteSt Plate « + «c+. « o 


C. Plate Motion Measurement. . 


D. Total Head Fluctuation Measurements. 


Lil. THEORY 


A. Forcing Moment Coefficient .... . 


Bo  Uumeilencesrramnes. « 2. .« « «© « e 


C. The Wave Analysis System . « « « « « 


IV PROCEDURE 


A. Total Head Turbulence Probe Calibration. 


B. .Tumnel Measurements. .« « « « «© « « 


Gis Data Analysis. e e @ + e ® e e e e 


V PRESENTATION AND DISCUSSION 


VI SUMMARY AND CONCLUSIONS 
VII RECOMMENDATIONS .... 
VIII REFERENCES. . 2. 2 ee 


APPENDICES 


ne Illustrations. ee 


B. Experimental Results 


C. Power Density Spectra. 


OF RESULTS. 


page 


33 
43 





18: 


i: 


LIST OF FIGURES 


Water tunnel test section showing electrometers, amplifiers, 


DC filters, and strain gage controls « « « « « « «© © 6 «© « e : 


Left: Sanborn 2000 tape recorder 
Center bottom: Precision Instrument Company portable tape 
recorder 


Right: Technical Products Company wave analyzer system. . « -» 
Proves afid supporting structuré. . . « «sx esseseveee 
Probes and supporting structure. . . » »+«eeeesese 0 @ 
Prove 7A" sigmal calibration . « «+» «eee s sess 8 8 @ 
Probe "B™ signal calibration . . .».+«-sseer eee ese s 


Probe: A" sigmal: velocity = O"Ttfeec. . 1s se ee ew ww ow 


Probe "B" signals velocity = 0 ft/sece « «ee we eee ececes 


Probe A‘ signal; velocity 1a ft/sec eee @ © © © © © © @ 


Probe Be signal; velocity 11 ft/sec “e377 ef © © © © @© 6 @ 


Bletesstenaieavelocity — Ouft ec, .)s «ss 6 oe 0% 2 6 


Platemsienal« velocity = 6G ftygec. « « « es + se ee ee ee 


Plate Signal; velocity Jl ft/sec eee e «© © @©@ © © © © © © 


Plate simnaly velocity — 14) ft/sec . «+ 6 + 0 6 ote 0 ee 8 


Turbulence probe arrangement, cross section of test plate and 


Coordimate SySiwme . 1. 2s se es we ew me wwe ewe ee tet lt 


Block diagram of instrumentation for one of the two turbulence 


comb en si eels Wes) 6 6 © © © 6 © #08 6 © 6 6 6 © 6 # @ © « 


Normalized cross power density vs. non-dimensional spanwise 


Clee meet eet eMecc Tels fe soos 6: 68-0 0 0 0 6s oe 0-0 
Correlation phase angle vs. non-dimensional spanwise distance. 


Normalized cross power density and vortex strength vs. self 


SxXCiicaiwon parameuer 6. « «5 « 6 4s :6) 6s #6 6s s « « « « 


page 


30 


30 
30 
30 
sil 
31 
31 
1. 
31 
ba 
32 
32 
32 
22 


eo 


wy 


39 
36 


aC 


203 


Zee 


Forcing moment coefficient vs. self excitation parameter. . « . 


Reproduced results; normalized cross power density and 


correlation phase angle vs. non-dimensional spanwise distance . 


page 
38 


SN) 





LIST OF TABLES 


page 
Tape Pecrstre dmonont UL@S« 6 « 6 6 6 6 © 6 «© 6 © seo. « os Al 


Table 2 [Pea etmONamtILIeCSs » © « 6 © « © 6 © 6 6 up e © 6 8 6 [Uomme L2 





I INTRODUCTION 


A. Flow-Induced Vibrations 

Structures subjected to a flow field may, under certain circumstances, 
experience flow-induced vibrations. These vibrations are undesirable and 
indeed may cause complete collapse of the structure due to fatigue induced 
by the vibrations. 

Unstreamlined or "bluff" bodies in a flow field shed discrete vortices 
in the wake. The shedding of these vortices gives rise to periodic lift 
forces in a direction transverse to the flow direction. It has been shown 
in reference (2) that for bodies of finite thickness this vortex-induced 
lift force increases linearly with the transverse vibrational velocity, 
thus producing what may be called a "self-excited" vibration. The fluid 
mechanism involved in this feedback between motion and forcing force is 
unknown but may be expected to be due to an increase in the strength of 
- the discrete vortices being shed and/or to an increase in the spanwise 


correlation of the vortices being shed. 


Be Scope of the Investigation 


The purpose of this investigation is to determine the cause of the 
amplitude dependence of the forcing moment so that engineers may devise 


steps to prevent the buildup of large amplitude vibrations. 





Re 


Il EXPERIMENTAL EQUIPMENT 


A. Water Tunnel 

The experiment reported herein was carried out in the closed-jet, open 
circuit water tunnel of the Hydrodynamics Laboratory at MIT. The water 
tunnel is fully described in reference (1). Figure 1 shows the water tunnel 
test section. Pump number III (rated capacity 3000 G.P.M. at a head of 45 
ft. and 1175 R.P.M.) was used. 

Flow measurement was obtained using the tunnel contraction as a flow 
meter. Its calibration was checked and found to be the same as reported in 


reference (2). 


B. Test Plate 

The test plate used was plate D from reference (2), having the 
geometry shown in Figure 15. The plate was mounted vertically in the test 
section in such a way as to permit rotation about the center of its semi- 


circular leading edge. The plate is fabricated of aluminun. 


C. Plate Motion Measurement 

1. Strain Gages 

Tatnall type C9-121-R2C strain gages were cemented to the torsion 
spring which restrains the plate motion and were waterproofed. The four 
gages form a wheatstone bridge which may be balanced by means of an external 
variable resistance. The bridge output voltage is proportional to the angle 
of plate rotation. See Appendix C of reference (2) for theory of calibra- 


tion of the strain gage bridge. 


2. AMplit ter 
A decade amplifier model 220-A made by Ballantine Laboratories Inc. was 


used to amplify the strain gage bridge output. 








3. Oscilloscope 


An Eico model 460 DC-wide band oscilloscope was used for visual 


observation of the strain gage bridge output. 


D. Total Head Fluctuation Measurements 

1. Total Head Turbulence Gages 

The gages for measuring the total head fluctuations (turbulence) are 
fully described in references (3 and 4). Figures 3 and 4, Appendix A, are 
photographs of the total head gages and Figure'15, Appendix A, is a drawing 


of the probes and test plate. 


2e Hlectrometer Preamplifier 


Two electrometers (one for each turbulence probe) were used to match the 
probe and amplifier impedences. 
The electrometers are fully described in references (3 and 4). 


Figure 16 is a block diagram showing all electronic components. 


3. Amplifiers 

Two Dana model 2000 solid state, wideband, differential, low-level D.C. 
amplifiers with integral power supplies were used to amplify the turbulence 
probe signals. The amplifier has a flat response to 73db from DC to 


approximately 70 KC. 


4. Oscilloscope 


A Dumont type 333 dual beam oscilloscope was used for visual observa- 


tions of the two turbulence probe signals. 


>« Portable Tape Recorder 


A Precision Instrument Co. portable tape recorder model PS 200 A was 


used to record the turbulence probe signals. 








Le 


6. Endless Loop Tape Recorder 


A Sanborn model 2000 endless loop tape recorder was used to give a 
continuous input of the turbulence probe signals into the wave analysis 


system. 


7. Wave Analysis System 
To analyze the turbulence probe signals, a Technical Products Co. 


TP 625 Wave Analyzer System was used. This system consists of the following: 


a. One TP 626 Oscillator 

This oscillator produces een frequency which is continuously 
variable from 97 KC to 122 KC over three different frequency bands: 
97.0 KC - 97.25 KC; 97.0 KC - 99.5 KC; and 97.0 KC — 122.0 KC. The 
oscillator output is used to modulate the signal to be analyzed in 


the TP 627 analyzers. 


be Two TP 627 Analyzers 

These analyzers determine the frequency and amplitude of the 
individual components of a complex wave within the frequency range 
of 2-25,000 cycles. The results are measured in percent of total 
signal. 

These analyzers are equipped with two matched TP-218 type 
filters»as»required for opeRauion in conjunction with a TP-645 


Multiplier. 


ce One TP-645 Multiplier 
This multiplier, when used in conjunction with the other components 
of the system, gives the cross power spectrum of two complex waves with- 


in the frequency range of 2-25,000 cycles. 





III THEORY 


A. Forcing Moment Coefficient 
In reference (2) Noutsopoulos develops the equation of motion of a 


single degree of freedom plate spring system vibrating about an axis through 
the leading edge of the plate. The motion is caused by a forcing moment 
resulting from flow separation at the trailing edge of the plate. Assuming 


two-dimensional motion, the equation of motion of the plate may be written 


I, + Co + Ko = M 
where I, = the total moment of inertia 
C, = the total damping coefficient 
K, = the total stiffness of the system 


iwt 
and & = the angular plate displacement in radians A= A, a : 


From purely dimensional considerations, the simplest form of the 


forcing moment can be written 
~~ = 22 
4 M 3 p H(2b) V 


where Mp is the spanwise average forcing moment coefficient, H is the span 
and V the ambient fluid velocity. When the forcing moment is periodic as 
in the case of excitation by vortex shedding, the moment coefficient is 


conveniently written in complex form as 


_ i(wt + @ ) 
=f 08 © 


where w = angular frequency of vortex shedding from one side of the plate 


in cps 


2 
I 


phase angle tetween the spanwise average forcing moment and the 
motion ef. 


t = time in seconds. 





6. 


The moment on the plate may also be written in terms of a lift force 
and an appropriate moment arm where the lift force is a function of the 


circulation which is a measure of the vortex strength. 


Therefore 
M ~~ 2bN 
where 2b is the chord length and 


N is the lift force. 
The lift force may be written as 
. 2 =p VA 
where J* is the circulation on the plate and is a function of position along 


the span and time such thats: 


H 
i = pv f ik (z,t)dz 
0 


If the circulation is considered to be a periodic function then 


["(2,t) = pe? gilt +0 + 8(2)) 


with 9 (z) being the phase angle of the circulation at some position, z, along 
the span, measured with respect to the value, ae at midspan of the plate so 


that 


H H 
N = pv oi (wt + @ ) J. [* o(z)eos Odz+i e Co (2) sin Qs 


Assuming application of the lift force at the trailing edge the 
moment on the plate may be written as 


mM == 20N 


H H 
M = 2p Vvbe"\™* bi > ) f | ((z)eos Odz+i [ Totersin 8 dz 
O 0 . 


Te 


It is now assumed that 


o~ f 


where Oe = circulation in the vortex 





us = the longitudinal fluctuating velocity component as measured by 
probe B 


and f = vortex shedding frequency. 


From reference (4) 


liar 


where i. = the velocity of the separation streamline. So that the assumption 


ro 


here is that 
? 
"3 w V. 


Substituting for hae the moment may be written 


: Hi caaese a 
M ~ pVbe" +@ ) f “ cos 9 dz + i { -B sin 8 dz 
QO 0 


Now assuming the moment coefficient to be a real quantity and equating 


the two expressions for moment we have 


= h H 
Re cs (rouv | J upp cos B dz 
‘ O 


and 
a al 
f ut= sinQ@dz = O 
0 B 


Here the assumption is that the wake is symmetrical behind the plate, so 
that the integral of an odd function is zero. Therefore, the forcing moment 


coefficient is proportional to the integral 





f uf* cos 9 dz 


where Th is the longitudinal fluctuating velocity component measured by 


the traversing turbulence probe (probe B). 


B. Turbulence Probes 

Reference (4) describes the development and theory of the total head 
tube in detail. Briefly the following is a resume of the theoretical con- 
siderations involved in the measurements. 

In general the instantaneous response of a total head tube may be 


written 


a = Ba Ply ture, (sve, ewe 
, ¥ 


22 2g 28 


where p is the pressure and u, v, and w are the velocity components in the 
coordinate direction. The barred quantities are mean values and the primed 
symbols are fluctuating values. 

From the geometry of the probe, it is a good assumption that the 
turbulence probe is much more sensitive to the velocity components parallel 
to the mean flow than to those which are perpendicular to the mean flow 
(i.e., the terms v'*/2g and w'*/2g are neglected). Also low frequency 
response of the piezoelectric crystal makes the turbulence probe insensitive 
to mean quantities. 

Therefore the instantaneous total head response becomes with these 
assumptions: 


_ pt re 2uut zt te 


tg g 28 28 


The D.C. portion of this signal (the time average of the fluctuation) 


is given by: 





9 


The D.C. level is removed by a filter so that 


Ho - Ho = Rhy wet, ule _ ule 
ee te s g 2828 


An important assumption at this stage is that all terms on the right 


hand side of the equation can be ignored with respect to uut/g, whereupon 


and the electrical signal as measured by the turbulence probe is proportional 
to the product of uu!, the product of the mean longitudinal wake velocity and 


the longitudinal fluctuating velocity component. 


C. The Wave Analysis System 

The two electrical signals measured by the two turbulence probes located 
in the wake of the plate and positioned along the span were analyzed by the 
wave analysis system. | 

In general the wave analysis system accepts one or two periodic or 
random signals as inputs and provides as output the power density of either 
Signal and the cross power density expressed as a complex number. The 
power density is a representation of the amount of power in the input signal 
contained in a given bandwidth. The cross power density is a measure of the 
power (contained in a frequency bandwidth) common to the two input signals. 

If the probe signals were harmonic they could be written as functions 
of time as f(t) and f(t) and their transforms as F, (nj) and F,(n2)- If 
F, (ny) is multiplied by Fa(n2), the complex conjugate of F,(n2), and the 
product is averaged in time, only products of common frequencies remain so 

(n) ———— (n) 

hie b ne F,(n)F,(n) = cross power density. ® 


AB ay also be 
written 
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® .(n) = V PsD, PSD 


B 


This can be shown since 


T 
F,(n) = z J fk? gin Wt st 
Ty 
F,(n) = . f fat) e7 inW,t a 
4 
and 
4 T, . 
: | Pate t, (+ tee at = i. { f, (@)at i +7) 
-Y 7 
1 
= ¥ pine ta 4 J fi (wee ol dt 
Nt—-cO 
yy 
therefore 
4 re) 
4 J fy (t)ea(t +F)dt = ba [Fine (n)| ginw,T 
Ty 1 2 : : 
4 Ut —0d 
and 


n % 
F,)(n)F2(n) = 4 : e inant at if f,(t)fo(t +7 )dt 
-T -T, 


but 





ll. 


% 
@ = 2 f e Nw? 4 5 J P(t it seat 
AB 1 aly } “Ty 


and since machine operation gives 


Fy(n) = PSD, 
F2(n) = V PSD, 


then 


d PSD i‘ V PSD, 


AB 
The TP 625 wave analyzer system accomplishes this by passing the random 
input signals £,(t) and f(t) through a set of amplitude and phase matched 
band pass filters before multiplication. In order to preserve the phase 
angle © between the two input signals, two modes of operation of the 
multiplier yield the real and imaginary parts of a complex number called 


the co-spectrum and the quad-spectrum respectively so that 


2, 


4 
| (co-spectrum)? + (quad-spectrum)?\ 


{ 


and 


Q i -l quad-spectrum 
= tan 
co-spectrum 


For random signals these quantities have the units of power density or power 
divided by the filter band-pass at the center frequency of the filter. 
Now if it is assumed that the signal from probe B is correlated 


(i.e., similar) with the signal from probe A, the fixed prove, we can write 


a = t 
ul K(z) ul 


where K(z) represents the proportionality between the signals as a function 





12. 


of their spanwise separation. As noticed previously the output of the 
i 11) %1)2 

probes are WU #0 that|, | ~ usugue. Note that it is tacitly 

assumed that u 


= Upe Therefore 








2 
~ utfuty2 . ututy2 
> | AB AB 


and if we divide by > |= we have 
AA 











2 
“) ufuftu2 ° ututu2 
AB AB 
Ree ee 
aC CHEE: 
Substituting for uy 
2 
“ K(z)ut? ° K(z)up 
aol 
2 
u f 
>. K(z UZ 
1te 
oNd UB 


Since wave analyzer operation gives 


b = PSD, 


AA 
then 
\$ |" 
ee NY ae 
— eae 
PSD, ue 


and referring to the discussion of the forcing moment coefficient where 





13. 
H 
-1 Te 
he w~w [ cov | { ug cos 9 dz 
0 


A 


are available as outputs of the wave analysis system and values of u are 


This suggests a method of evaluating the integral. \d and PSD, values 
AB 


easily measured in the tunnel. Values of @ are available as shown from 


Q ei ae quad-spectrum 


co-spectrum 





lhe 
IV PROCEDURE 


A. Total Head Turbulence Probe Calibration 

The total head turbulence probes, electrometer preamplifiers, D.C. 
filters, and the turbulence probe signal amplifiers were calibrated as 
assembled units. The two system outputs were adjusted by changing the 
amplifier gain on one system so that the output of the two systems were 
equal for identical inputs. Figures 5 and 6, Appendix A, show the 
amplified turbulence probe signals when the probes were subjected to a 


- step impulse. See reference (2) for calibration details. 


B. Tunnel Measurements 
The fixed turbulence probe was initially positioned in the tunnel at 
y/t =O and x/t = 1. See Figure 15, Appendix A for description of coordinate 
system. Alignment of the traversing turbulence probe (probe B) with respect 
to the fixed probe (probe A) is assured by the method of attachment of the 
_ two probes and the proper positioning of the traversing probe seal at the 
‘bottom of the test section. Measurements were taken with the plate free to 
vibrate at ambient velocities of 8, ll, and 14 fps, and clamped in a fixed 
position at an ambient velocity of 14 fps. The velocities were obtained 
from the contraction meter as described in reference (2). Mean wake 
velocities were measured as described in reference(3). While maintaining 
a constant ambient velocity by operation of the inlet and foot valves, the 
turbulence probe signals, plate vibration signal and a marker signal were 
recorded on four channels of the protable tape recorder for each spanwise 
position (z/t) of the traversing probe. Prior to and during the actual 
recording, the signals were monitored on oscilloscopes to insure that the 


proper signals were recorded. The pressure in the tunnel was maintained 
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at the same value by operation of the foot valve and was kept sufficiently 


high (3 psi) to prevent cavitation. 


C. Data Analysis 

The recorded signals were analyzed with the Technical Products Company 
TP 625 Wave Analyzer System. Each run consisted of a recording of the four 
signals (two fluctuating total heads, plate vibrational amplitude, marker) 
at a particular z/t and velocity. For each run the probe signals were 
rerecorded on a Sanborn-Ampex 2000 continuous loop tape recorder for con- 
tinuous play back into the wave analyzer system. Prior to each analyzation, 
the wave analyzer system was carefully calibrated with a 100 cps, 0.46 volt 
peak to peak peaiiator Signal. The gains of the two wave analyzers were 
set the same and the oscillator (TP 626) output indicator was adjusted with 
a reference oscillator to ensure accurate frequency readings. The gains 
of the two analyzers TP 627 were matched and the Power Spectral Density set- 
tings of the multiplier TP 645 were adjusted to give the same output for 
identical inputs. In addition the phasing between the identical input 
signals was zeroed using the phase zero dial. Calibration in this manner 
ensures that the relative power levels in the signals are preserved in the 
output, and that zero phase angle is indicated for identical input signals. 
The output of the wave analysis system is in volts and is only a relative 
measure. A 50 cps. band width filter and a long averaging time (2 seconds) 
were chosen for the analyzation. The filter and 2 second averaging time 
were selected to give smooth curves from which data points could be 
accurately determined. 

For z/t = 1 and 9 at each velocity, the power spectrum density for 
each signal and the co-spectrum and quad-spectrum were obtained. For the 


remaining valves of z/t, the power spectrum density of "B" (the movable 
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turbulence probe), co-spectrum and quad-spectrum were obtained. 

The horizontal scale of the x-y plotter was adjusted so that a full 
sweep covered from 50 to 200 cps. The vertical scale of the x-y plotter 
was set to 1 volt per inch. Also for z/t = 9 at each velocity, the power 
spectrum density "B" was obtained using a 1 cps band width filter and an 
averaging time of 0.2 seconds to indicate the exact frequency of the 
energy in the wake. Similarly for the plate motion signal, a spectrum was 


obtained to pinpoint the frequency of plate vibration. 





17. 
V PRESENTATION AND DISCUSSION OF RESULTS 


The outputs of the Wave Analyzer system are the power density spectrum 
of a single signal, the co-spectrum, and the quad-spectrum of the two signals 
taken jointly. These are presented for each run as Appendix C. 

For analyzation purposes, frequencies of 75 cps for runs with the test 
plate free to vibrate and 105 cps for runs with the test plate fixed were 
chosen. These frequencies were selected because they were the respective 
points of peak power as indicated by the plots of power spectrum density 
(PSD) and hence represent the frequency of vortex shedding. This reason 
for frequency selection is reinforced by the PSD of the vibrating plate 
Signals which also has its peak power at 75 cps. See the power spectrum 
density for vibrating plate signals, runs ll, 16, and 21, Appendix C. 

Power density spectra were obtained from the wave analyzer system with 
narrow filter settings as well as with wide filter settings for runs 6, ll, 
16, and 21. These narrow filter power density spectra show more clearly 
the exact frequency of peak power than do the power density spectra made 
with wide filter settings. 

Table 1, Appendix B, lists the values of PSD, co-spectra and quad- 
spectra used in the following analysis. 

Power density spectra for probe A (fixed probe) are shown in Appendix C 
for runs l, 6, 7, 11, 12, 16, 17, and 21. These are all runs with probe 
separation (z/t)= 1 and 9. It can be seen that for each velocity PSD, is 
the same for all probe separations z/t. This indicates that at small probe 
separation, one probe does not affect the signal of the other probe and 
that the signal from the turbulence probe A is stationary. 

The power density for the fixed probe A increases in magnitude at the 


plate vibration frequency with an increase in ambient velocity. However 
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the power densities for the fixed plate (runs 1 through 6 of Appendix C, 
ambient velocity - 14 ft/sec) are higher than are the power densities for 
the vibrating plate runs with the same ambient velocity as the fixed plate 
runs (see runs 17 through 21 of Appendix C for vibrating plate runs with an 
ambient velocity of 14 ft/sec). Thus power density is not a function of 
ambient velocity alone. 

If one computes the self excitation parameter of reference (2), 

.- = ee 5 rer each run made 5 this study and compares om and V*, 
it can be seen that Ors decreases as the self excitation parameter V* 
increases. V* is the ratio of the transverse velocity of the trailing 
edge of the test plate to the ambient water velocity. 

It has been shown in Section III that the power density spectrun, AA? 
is proportional to (utu)®. Thus ®, . / is proportional to (ul )?. The 
quantity, (u')*/f, may be considered a measure of the strength of the 
individual vorticies. Thus b 


A 


vortices being shed. Figure 19 shows the strength of the vortices 


et is proportional to the strength of 


® fer being shed plotted vs. the self excitation parameter V*. Figure 19 
shows that the vortex strength drops off rapidly with an increase in the 
self excitation parameter V*. Thus the increase in forcing moment co- 


efficient Mh. with an increase in the self excitation parameter V* cannot 
fe) 


be explained by an increase in the strength of the vortices being shed. 
For the case of the fixed plate (V* = 0) it appears that the values 
of vortex strength and coherency may vary with the ambient velocity. 
The power density spectra for the signal from probe B at small 
separations (runs l, 7, 12, and 17 of Appendix C) are essentially the same 


as the power density spectra for the signals from probe A. 





OF 
As shown in Section III: 


\§ aB\ = | (co-spectrum)? 1 (quad-spectrun)?| 


Op - igere quad-spectrum 
co-spectrum 


” 


Table 1 Jists the values of | asl and § for all runs. As explained in 
Section IV, Procedure, the wave analyzer system was adjusted so that the 
outputs were all relative to the 100 cycle sine wave calibration signal. 

The Technical Products Company wave analyzer system when in the power 
spectrum density mode of operation, multiplies the input signal by itself 
and in effect becomes a squarer. Time averaging then eliminates all 


products except those of common frequencies. Thus the power spectrum 





density of signal A = Om ~ (uta)? likewise @ 25 ~ (utu)?. When the 
wave analyzer system is used to obtain the cross power spectrum density 
the two input signals are again multiplied and time averaged. 
When interpreting the machine output we must remember that all d 's 
-are averages over a finite bandwidth. Thus when we look at a or Pop 
with the filter centered at frequency no aa) or F,(n,) may be zero yet 
>. or ®.. can be non zero due to contributions at other frequencies 
within the filter pass-band. When looking at d.. we are most interested 
in its value at the vortex shedding frequency. If this frequency is only 
slightly different at the locations of probes A and B, and the random 
turbulent background at surrounding frequencies is negligibly small, © 5 
can be zero and yet Ban om and Qo. would be non-zero! 
= a3 mn (atuta ie oo). Das ~ (ua) u)? and O33 ~ ~ (uta), then 
\O as \~ (Os aa| “16 BB Mere us and ug have the same frequencies. 
Therefore \$ aBl is a measure of overall signal strength as well as a 


measure of the similarity of the two signal amplitudes. To eliminate the 
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effect of overall signal strength all values of | in for each velocity 


can be normalized by their respective ® values (for each velocity 


~ AA AA 


is stationary as shown in Appendix C). 

Therefore: a = \@ | / d aq and is only a measure of the 
similarity of signal A and signal B. Values of }™ are tabulated in 
Table 1, Appendix B. Figure 17 is a plot of o* vs. z/t for different 
values of the dimensionless self excitation parameter V* = 20 UNche, ° 
Figure 18 is a plot of @ vs. z2/t for different values of V*. 

Vibrating plate runs at an ambient velocity of 14 fps were repeated 
after a two week interval to investigate the reproducability of the 
experimental data. Figure 21 compares values of b* and OQ vs. Z/te de 
can be seen from this figure that the cross power spectra compare very 
favorably, but a phase angle variation is to be expected. 

From Figure 17 it can be seen that the correlation of the turbulence 
probe signals for vibrating plate runs decreases as probe separation 
increases which is what one would expect. Figure 17 also shows that the 
normalized cross power density increases for vibrating plate runs with an 
increase in the self excitation parameter V*. 

Figure 19 shows that the normalized cross power density Q” increases 
as the self excitation parameter V* increases. 

Figure 18 shows that for z2/t< 5 for vibrating plate runs there is a 
gradual increase in © as z/t increases, but there is no significant change 
261 e for an increase in the self excitation parameter V*¥, For probe 
separation z/t > 5 for vibrating plate runs @ increases with an increase 
in z/t and decreases as the self excitation parameter V* increases. 

Figures 17, 18, and 19 show that the shedding of the vortices bec eneS 


more correlated and more in phase as the self excitation parameter V# 


increaseSe 


a ie 
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For the self excitation parameter V* = O (fixed plate) with a velocity 
of 14 ft/sec, Figure 17 shows >” to decrease to near zero at z/t w7 and 
b* to increase for values of z/t > 7. For V* = O Figure 18 shows 9 to 
increase gradually for values of z/t < 6 and 9 to increase rapidly to 
values > 180° for values of z/t 4 6. This indicates that for V* = 0 
there may be a spanwise cellular structure with alternate cells 180° out 
of phase. Assuming symmetry the wavelength of the cell would be 14 body 
thicknesses. Roshko (6) reports a 18 diameter wavelength for measurements 
behind a circular cylinder at a body Reynolds number of 80. 


As shown in Section III: 


cos 9 dz 


H | \* 
n. pies (DHEVa,? )~* f \@ asl 
i 0 BD as 
Values of (\d in cos 9 )/( ® An ) were plotted vs. z/t and inte- 
grated from z/t = 1 to z2/t = 9. It was assumed that the value of this 
integral is proportional to the integral evaluated over the entire plate 
Spane 


was calculated by 


o 
2 2 
Me ow (eve? f Baal? cos d(z/t) 
0 1 Qaa 


These values of Mp for the vibrating plate runs are plotted vs. V# in 
fe) 


Figure 20, Appendix B. The dotted line in Figure 20 is a plot of m 


af 


f VSe 


fe) 
V* as determined from reference (2). It can be seen that these two plots 


of m, are in qualitative agreement. In an effort to determine any pre~ 
te) — 
dominence of either lS aal?/ ?.. or Q in the integrand of the equation 


ie 
fe) 


herence, the predominant factor in 


used to determine m, (i-e., is the amplitude coherence or the phase co- 





22 


2 
De Ww (rvuz)7+ | \3,{ cos 9 d(z/t)? ) it was assumed that cos 9 was 
Oo 4 ] AA 


constant, hence 


z/t=9 2 
mow ceviey? Gad d(z/t) 
= if ® an | 


The resulting Tp plotted vs. V* is shown in Figure 20. 


2 
Next it was assumed that {6 | )/( ®. ) was constant, hence 


79 
.™ (fVu2) |cos 9 d(z/t) 
O 
.t 


m 


The resulting m, plotted vs. V* is shown in Figure 20. 


f 
fo) 
The conclusion drawn from this fugure is that it is impossible to 


determine the predominate factor in 
2 
Me ow (rvuz)7? fae cos 9 a(z/t) 
es ® AA : 


- 2 
Undoubtedly both (\ : )/( do) and cos © are importante Indeed it may 


A 
be said that amplitude correlation and phase correlation have no significance 
unless presented together. Phase correlation has no significance unless 
there is some amplitude correlation and vice versa. 

The- calculated values of Mp for the fixed plate runs assuming that 
the integrals evaluated from yee 1 to z/t = 9 were proportional to the 
same integrals evaluated Bret the entire plate were not considered valid. 


However, by extrapolation of the integrand over the entire length, assuming 


a cellular structure with alternate cells 180° out of phase, values of Le 
O 


obtained are in qualitative agreement with reference (2). 
The co-spectrum of the vibrating plate runs, runs 7 - 21, Appendix C, 


show that some correlation exists at 140 cycles. This correlation was not 
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investigated, but it could be the correlation of the first harmonic of the 
vortex shedding frequency. The vibrating plate spectra for runs lb, 21 
Appendix C, show power at about 140 cycles, and this could explain the 


correlation at 1,0 cycles. 
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VI SUMMARY AND CONCLUSIONS 


Spanwise cross correlations of the longitudinal turbulent velocity 
component in the early wake of a flat plate were measured in water under 
conditions of flow-induced torsional vibration of the plate in order to 
ascertain the nature of the known dependence of the forcing moment upon 
the ratio of vibratory plate velocity to free-stream fluid velocity. 

The following conclusions are drawn: 


1. The increase in the forcing moment coefficient Mp with the 
Oo 


2 self excitation parameter V* is due to an improvement in correlation 
or coherency of the wake along the span rather than to increasing 


strength of the individual vortices being shed. 


2- The forcing moment coefficient can be expressed in the form 


"Saal 
Me a (fbHvu2 7? | 2 sal cos 8 dz 
o Pay 
; + 
3- Both the absolute magnitude of the cross power density ® 
and its phase angle Q are important in explaining the relationship 


between the forcing moment coefficient m, and the self excitation 
fe) 


parameter. Spanwise variations in phase angle are shown to be 
reduced with increasing V* and coherence of vortex strength at the 


vibrational frequency is increased. 


4. Measurements of the spanwise wake structure under conditions 
in which plate vibrations are prevented are consistent with a cellular 


structure with a phase shift of 180° between alternate cells. 


5- In order to minimize self-induced vibration, engineers should 


design structures in such a way as to eliminate wake coherency. This 


256 


might be accomplished by a random spanwise variation of trailing edge 
geometry, or by attaching to the trailing edge short independent 


sections which are free to rotate about the trailing edge. 
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VII RECOMMENDATIONS 


1. In the co-spectrum of probe B some indication of correlation was 


observed at 140 cps. It is recommended that this be investigated. 


2 It is recommended that measurements such as taken in this study 
be taken so as to include observations of the wake along the entire span. 


This will necessitate a new turbulence probe supporting mechanism. 


3. It is recommended that a 5, 10, or 20 cycle bandwidth filter be 
obtained for the Technical Products Wave Analyzer System. These more 
narrow filters will permit a more accurate determination of the frequencies 
of peak power in both the power spectrums and the cross power spectrums. 
The available one cycle filter proved to be too narrow for the random 


Signal involved. 
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TITLES OF PICTURES ON THE FOLLOWING PAGE 


Figure l: 


Figure 3: 


Figure 43 


Water tunnel test section showing electrometers, 


amplifiers, DC filters, and strain gage controls. 


Left: Sanborn 2000 tape recorder. 


Center bottom: Precision Instrument Company portable 


tape recorder. 


Right: Technical Products wave analyzer system. 
Probes and supporting structure. 


Probes and supporting structure. 








SATLIL 
4unLorld 
YO 
HOVd 
SsnorAddd 
das 








o 


—- > 
—_- 
- aS 
== 
= eS a Se 
= ? 
— i 
—_ 
: a T 





PGs FIG, 6 
PROBE A SIGNAL PROBE B SIGNAL 
CALIBRATION 


CALIBRATION 





ees | FIG. 8 
PROBE A SIGNAL PROBE B SIGNAL 
VEL.= O FT. /SEC. 


VEL.= 0 FT./SEC, 





FIG. 9 FIG. 10 
PROBE A SIGNAL PROBE B SIGNAL 
VEL.= 11 FT./SEC. 


VEL.= 11 FT./SEC. 








r ar earenr 
A. PANAMA. Ws 


ee ea 
ani SEi.2ee 
1aaee Pee 
aT 





FIG. 11 FIG. 12 
PLATE SIGNAL PIATE SIGNAL 
VEL.= 0 FT./SEC. VEL.= 8 FT./SEC. 





FIG. 13 FIG. 4 
PLATE SIGNAL PLATE SIGNAL 
(ass Thy Suey Sale VEL.* 14 FT./SEC. 





~~ 2 oo ee ee 


sis ~ é _ i > Li SE Ri nt NT eet Ow 4 Oe He oo 





-— as 


Z 


The origin of the coordinate system is located at the 


center of the test section at the trailing edge of the test 


plate. 


Figure 15: 


Turbulence probe arrangement, cross section 


of test plate, and coordinate system 








———— TURBULENCE PROBE 


OSCILLOSCOPE 


PREAMPLI FIER 


FIG.16 


arb iin 


A PLIFIER 


FLOCK DIAGRAM OF INSTAUMENTATION FOR ONE OF THE TdO 


TURBULENCE PROSES 





3h. 








APPENDIX B 


EXPERIMENTAL RESULTS 








356 


aoueystp estmueds 


TeuotsuseuTp—-uou *sa AYTsuep Zemod ssozo poeZzTTeULIoy 


+t 
z 
6 va 5 c | 
CSOO=« A 0-0 
ILOO=¥ AX 
LLOO=2 A fe’ Q 
O Yy 
© 
: 7 
Oo 
O=2A e 


JINVLISIG 3SIMNVdS= 
YSLAIWVYVd NOliVLIDX3 3135S 


“3 


NOILVI3dNu09 SSO) 40 Se NOE or 











2/4 eansty 


20 





ee Oe ee eee ES ee a eae ee ee ey SS ee 


240 





@=PHASE ANGLE OF CROSS CORRELATION 
V*=SELF EXCITATION PARAMETER 


te f =SPANWISE DISTANCE 


‘> V* =O 
200 


a / 


I60 
Ov *=0052 
140 


120 


100 


C—O ees NNT 


80 


60 


0 AV *=007I 
& 


a V *=0077 


40 ue 
ee 


© 
e CO} ee 
20 


x 
4 

a 
9 


\ 





chy 


Figure 18: Correlation phase angle vs. non-dimensional spanwise 


distance 
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Figure 19: Normalized cross power density and vortex 


strength vs. self excitation parameter 
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Figure 20: Forcing moment coefficient vs. self excitation 


parameter 
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Figure 21: 


G° = MAGNITUDE OF CROSS CORRELATION 





and correlation phase angle vs. non-dimensional 


spanwise distance 


Reproduced results; normalized cross power density 
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Table 1 Measured Quantities 


Al. 





Run dimensionless 


Ou FWD BP 


z/t 


distance 


PSD 


(volts)* 


3025 


Lh 


La | 


coO=- 


spectrum 


1.55 
1.50 
aE) 
O45 


-0.55 
-1.25 


0.31 
0.22 
0.17 
0.10 
0.05 


0.90 
0.58 
0.39 
0.21 
0.05 


0.85 
0.53 
0.18 
0.05 


-0.10 


quad- 
spectrum 
(volts)* (volts)# 


0.12 
0.20 
0.33 


-0.10 
-0.17 
-0.73 


0.07 
0.06 
0.05 
0.05 
0.05 


0.15 
0.30 
0.25 
0.20 
0.20 


0.00 
0.21 
0.24 
0.17 
0.08 


oO x 10° 


(radians) 


plate 
fixed 


raat 


u 


(ft/sec) 


10.85 


ft 


V 
(ft/sec) 


$f 





# Values measured at 105 cps for fixed plate and 75 cps for vibrating plate. 





Table 2 Derived Quantities 
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Run — ae d AB 0 a ay x10" 18, .\cos 8 (o.,\ 
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at plate 1.55 pele 0.48 L007 6.260 0.0198 

2 — O 1.52 is 0.47 3.98 6.000 n 

3 1.19 16.0 0.37 3.13 3.560 " 

hy 0.46 67.5 Os 1.19 ~0.539 " 

5 0.58 198.8 0.18 Be ~0.832 " 

6 dusts 210.3 0.45 2581 Seezirie " 

7 0.077 Ows2 ey, 0.58 15.90 4.950 0.0150 

g " 0.23 1563 6.12 11.50 2-570 " 

9 " 0.18 16.4 Or52 8.80 1.550 " 
10 " 0.12 26.0 0.20 5.50 0.546 t" 
121 ‘" 0.07 45.0 Onis 3.57 0.173 " 
12 0.071 0.91 9.5 0.56 6.23 5.560 0.0181 
gy " 0.65 Piped 0.40 bol5 2.570 " 
A, " 0.46 32.7 0.28 sical 1-226 " 
15 " 0.29 1.3.6 0.18 2.00 Delay n 
16 " 0.19 Pac oneal ney hit 0.077 i 
19; 0.052 0.83 3.6 0.42 5.83 4.780 0.0278 
18 " 0.57 oy 0.29 nO 2.100 " 
19 " 0.30 53.0 0.15 2.09 0.375 " 
20 " 0.18 73.6 0.09 1.25 0.063 " 
21 " ov 146.3 0.08 pel 20.113 " 
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APPENDIX C 


POWER DENSITY SPECTRA 
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